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Abstract 



Current-voltage (/ — V) characteristics and the increase of superconducting 
£> ■ transition temperature T c as a function of the number of half-unit-cell layers 

n are studied on Bi2Sr2CaCu20s+5 thin films grown by molecular beam epi- 
taxy. These observations are interpreted in terms of Kosterlitz-Thouless (KT) 
ON ■ transition associated with vortex string threading the film. Considering the 



change in the types of the lowest free-energy excitations as a function of the 



a 

*^ • size of vortex pair, we have formulated KT renormalization-group (RG) equa- 

tions. They are with respect to independent vortex and antivortex pairs with 



distance r less than n-dependent characteristic length r n and then in terms 
of vortex-string and antivortex-string pairs with r > r n . RG equations of 
these two regions are connected by the continuity of superconducting carrier 
density and vortex density per area. Using experimental data obtained with 
Bi2Sr2CaCu208+(5 thin films, we have solved the RG equations. It is shown 
that the KT transition occurs in the vortex-string region with interstring lnr 
interaction only. Our calculation successfully explains KT scaling in I — V 
characteristics and the relation between T c and n. These results indicate that 
the interlayer interaction plays an important role in Bi2Sr2CaCu20§+5 thin 
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films. 
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I. INTRODUCTION 



Two-dimensional (2D) Kosterlitz-Thouless (KT) transition has been proposed as 
a possible mechanism for high-T c superconductivity. Anisotropically Josephson coupled 3D 
XY model served as a microscopic model for layered superconductors with weak interlayer 
coupling. In the anisotropic 3D XY model with infinite volume, Hikami and Tsuneto 
(HT) have proposed two types of excitations, i.e., independent vortex pairs in each layer 
and closed vortex rings. || This view is also supported by Korshunov's theory M based 
on the 3D Josephson Junction array. Superconducting transition in thin metal films, |7J 
Bi 2 Sr 2 CaCu 2 8+(5 (BSCCO) single crystals, f8HH YBaaCusOr^ (YBCO) single crystal, 
|12|1 and Tl 2 Ba 2 CaCu 2 0s (TBCCO) single crystal |13[ were interpreted in terms of 2D KT 



theory. 2D KT transitions in monolayered or thin-film YBCO have also been observed. 
fl4Hl6| In bulk samples, superconducting transition have been considered to occur owing to 
KT transitions by independent vortex in each Cu0 2 layer despite the presence of interlayer 
coupling. []8| p~3| 



In this paper, we will show by using thin films that excitations of vortex string threading 
the film do emerge as a result of interlayer interaction. Extending HT theory to thin-film 
cases, we are naturally led to the concept of pairs of vortex string and antivortex string 
threading the film. This concept has been introduced and demonstrated experimentally by 
Matsuda et al. using YBCO thin films prepared by reactive evaporation technique. |17| 
They observed the increase in superconducting transition temperature T c as a function of 
the number of unit-cell layers and argued that this increase in T c is due to vortex string and 
antivortex string. However, they did not consider the existence of independent vortex and 
antivortex pairs in their renormalizat ion-group (RG) analysis. This is not always justified. 

The purpose of this paper is to show that KT transition by vortex string occurs in BSCCO 
thin films even when we consider the coexistence of vortex-string and antivortex-string 
pairs and independent vortex and antivortex pairs. It is generally thought that BSCCO is 
more anisotropic than YBCO. Then, the role of independent vortex and antivortex pairs 



3 



seems more important in BSCCO than in YBCO even if we consider the KT transition by 
vortex string. Consequently, a more rigorous treatment considering the coexistence of both 
excitations seems necessary. Thus, to elucidate the role of interlayer interaction, we study 
the superconducting properties of BSCCO thin films grown by molecular beam epitaxy 
(MBE) and obtain the results to be interpreted by KT transition by vortex string. We 
observe the universal jump in the temperature dependence of nonlinear power law a(T) 
obtained from current- volt age (/ — V) characteristics and the increase in T c as a function 
of the number of half-unit-cell layers n for our BSCCO thin films. The increase in T c as 
a function of n is rather similar to that of YBCO thin films and can be regarded as the 
evidence of the presence of vortex strings threading the film. This is in contradiction to the 
expectation that for very anisotropic BSCCO film the interlayer coupling is so weak that 
independent vortices may be dominant excitations. We further confirm that the universal 
jump in a-T relation occurs in the vortex-string form and obtain the size of pairs of vortex 
string and antivortex string. This is analyzed on the basis of the KT RG flow with interlayer 
coupling. To consider the strong anisotropy of BSCCO superconductors, we formulate 2D 
KT theory for multilayered thin film coupled by interlayer Josephson coupling. This theory 
includes three types of low-lying excitations, i.e., independent vortex and antivortex pairs 
in each layer, vortex-string and antivortex-string pairs threading the film, and closed vortex 
rings as proposed by Hikami and Tsuneto. || Comparing the free energies of these three 
types of excitations, we find that the lowest free-energy excitation changes its form from 
the independent vortex into the vortex string as a function of r, the separation between 
independent vortex (or vortex-string) and antivortex (or antivortex-string) pair. We consider 
the two regions separated by a characteristic length r n . r n will be defined later in Eq. 
Then, the lowest free-energy excitations are independent vortex and antivortex pairs with 
distance r < r n , and are vortex-string and antivortex-string pairs threading the film for 
r > r n . Thus, we write KT RG equations with respect to the independent vortex pairs 
with In r intralayer and r 2 -dependent interlayer interaction for r < r n and then with respect 
to vortex-string pairs with lnr interstring interaction for r < r n . Differential equations 
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of these two regions are connected by the continuities of superconducting carrier density 
and the vortex density per area. In this way, we have reduced the 3D XY problem to the 
2D XY problem with interlayer coupling. Solving the 2D KT RG equations in BSCCO 
thin-film case, we obtain RG flow in which interlayer coupling is considered. We observe 
the transition from independent vortex to vortex string in RG flow at any temperatures. 
We shall see that KT transitions occur in the vortex string region in these RG flows thus 
obtained. Considering the asymptotic behavior of the free energy of a vortex pair we show 
that the r 2 dependent interlayer coupling necessarily leads to the KT transition in the form 
of vortex string threading the film. Analyzing the experimental results on both I — V and 
T c -n relations consistently, we conclude that the interlayer coupling plays a significant role 
even in the very anisotropic BSCCO so that the KT transition occurs in the vortex-string 
excitation. Considering the coexistence of the independent vortex, we show that the KT 
transition by vortex string occurs in films as thick as 20 half-unit-cell layer for the strongly 
anisotropic superconductors BSCCO. Comparing the results of RG flow with a-T relation 
we can determine the size of pairs of vortex string and antivortex string during the KT 
transition. Thus, we believe that our experimental and theoretical method of analysis using 
thin films will pave the way to a quantitative understanding of the interlayer coupling in 
superconductors such as BSCCO, YBCO, etc. This will be a step forward to revealing the 
nature of superconductivity in these materials. 

This paper is organized as follows. In Sec. 2, RG analysis of KT transition in thin films 
is formulated considering interlayer Josephson coupling. Using the parameters of BSCCO 
thin films we will show that KT transition occurs with respect to vortex-string pairs. In 
Sec. 3, the process of sample preparation and the method of measurements are described. 
In Sec. 4, the experimental results on I — V characteristics and the T c vs n relation are 
described and interpreted in terms of KT theory in which vortex string plays the role of 
vortex in pure 2D KT theory. In Sec. 5, the discussion and the conclusion are given. 
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II. RENORMALIZATION GROUP ANALYSIS OF KT TRANSITION IN THIN 

FILMS 



KT transition in 2D superconductors is related with binding and unbinding of quantized 
vortex and antivortex due to the logarithmic interaction between them. This transition oc- 
curs at a temperature lower than Ginzburg-Landau (GL) mean-field transition temperature 
T c0 , below which condensate wave function is formed. In high-T c superconductors, the effect 
of interlayer coupling should be considered. 

Now we consider three types of low-energy excitations in 2D XY model with Josephson 
interplane coupling. On the basis of anisotropic 3D XY model with infinite volume, the 
effect of interplane coupling is considered by Hikami and Tsuneto. f| According to them, 
the energy of an independent vortex and antivortex pair in each layer is given by || 

U vp (r) = <2>kK x ln(r/e„ 6 (T)) + 2vrK ± (r/rf) 2 + 2E C . (1) 

Here, K\ is the coupling constant in the ab plane, K± is the coupling constant along the 
c axis, £,ab{T) is GL coherence length in the ab plane at temperature T and is defined as 
£,ab(T) = £ a &(0)(l — T/Tco) -1 / 2 , d is the interlayer distance, and E c is the core energy to 
create a vortex in one layer. The second term in Eq. ([[]) expresses the interlayer Josephson 
coupling energy. In the presence of interlayer coupling, Hikami and Tsuneto showed that 
there are closed vortex rings within a sample. The energy of a vortex ring in an infinite 
system is given by || 

U ring {r) = (r/d^TcK, In(r/^(T)) + 2E C ). (2) 

r/d is, roughly speaking, the number of the layers which the vortex ring crosses. U vp (r) = 
U r i n g(r) yields a characteristic length defined as r = d(Ki/K±) ln(ro/£ &(T')). || Here, we ig- 
nore E c since it is small. The anisotropy ratio K\j K± is defined as K%/K± = (£ a b(T)/£ c (T)) 2 , 
where £ C (T) = £ c (0)(l — T/T c0 ) -1 / 2 . Then, (a) for r < r , independent vortex pairs in 
each layer dominate, (b) For r > r , closed vortex rings in a 3D volume dominate. 
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Application of HT theory to thin film leads to the concept of vortex string as follows. 
For BSCCO, we have r (T=64.5 K)=l xlO 4 nm, using the physical parameters in Table 
|. Then for our n-layered films (n < 20), the relation r$ > nd is satisfied. In this case, 
closed vortex rings cannot exist within the films. Then, there can be only pairs of vortex 
and antivortex strings threading the film and those of independent vortex and antivortex. 
The interaction energy of a vortex-string and antivortex-string pair in n-layered film is given 

by m 

U st (r) = n(2TtK x ln(r/U(T)) + 2E e ), (3) 

which is free from r 2 -dependent interlayer coupling term and less than U r i ng (r) for r > nd. 
Considering the entropy change for multilayered system, we have the free energy of the 
independent vortex pair as 

F vp (r) = 2vK x \n(r/U(T)) + 2irK ± (r/d) 2 

+2E c -k B T\n(n). (4) 

The fourth term is the entropy term for multilayered system which is to be understood as 
the difference from that of vortex string. The origin of the entropy term of Eq. @) will be 
explained later. Equating F vp {r) with U st (r), we introduce a characteristic length r n defined 
by 



ds 



n - l)K x ^ / r n \ + k B T\n(n) ^ 



\ K± \U(T)J K 
This is to be distinguished from r n defined from the relation U vp (r n ) = U st {r n ) by Matsuda 



et.al. |17j In Eq. fl5|), E c is ignored again. For r < r n the favored excitations are independent 
vortex and antivortex pairs, and for r > r n they are vortex-string and antivortex-string pairs. 
We also introduce another parameter l n = ln(r n /£ a &(T)). The values of l n (r n ) are listed in 
Table |. These values of l n are obtained by using renormalized K\ as will be described later. 

The transition from independent vortex pairs to vortex-string pairs can be formulated 
in 2D KT theory with interlayer coupling. We find the lowest free-energy excitation among 



the three is independent vortex pair for r < r n and vortex-string pair for r > r n under the 
following conditions, 

K 1 /K ± ^>K P , (nd<U(T)) (6) 

K 1 /K ± ^>K r , (nd>U(T)) (7) 

where 

K p = (5A36/(n-l))(U(T)/d) 2 , (8) 

K r = n 2 /((n - 1) \n(nd/U(T))). (9) 

For BSCCO film with n=20, we have K P =1.U, K r =9.U and KJK^IO 3 using the pa- 
rameters in Table Q. Then, the inequality K\/K±_ 3> K p and K r is satisfied. On this basis, 
we formulate 2D KT theory with respect to independent vortex pair with energy U vp (r) of 
Eq. ([]]) for r < r n and that with respect to vortex-string pair with energy U s t(r) of Eq. (|3]) 
for r > r„. Then, the differential equations of the two regions are connected. We will find 
that the KT transition occurs in the vortex-string region with interstring lnr interaction 
of Ugt(r). After solving the RG equations, we will reconfirm in Sec. 5 that the lowest free- 
energy excitations are such as we have supposed initially. A general proof for the absence 
of dissociation of independent vortex and antivortex pairs will be also given in Sec. 5. 

Now, we introduce in advance the various transition temperatures to be used in this 
paper. T^- T is the KT-transition temperature of a single half-unit-cell layer system taken 
as a scaling parameter for T. T c (n) is the superconducting transition temperature for n- 
layered system at which the measured resistance vanishes. T'Jji) is the superconducting 
transition temperature for n-layered system calculated from the RG equation. T c o(n) is the 
Ginzburg-Landau (GL) mean-field transition temperature for n-layered system. T m ^(n) is 
the midpoint temperature for n-layered system. T mid {n) is regarded as T c0 (n). These values 
are listed in Table |. 

The starting equations for our KT RG scheme are written as 

e(r + dr) = e(r) + 4n 27crdrP(r), (10) 

2kg 1 
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P(r + dr)=P(r)exp(- d p£-), (11) 



k B T / 

where dU (r) is given by 

dU(r) = dU vp (r) = + dr, 

(independent vortex) (12) 

dU(r) = dU st (r) = (-^-) dr ' 
\e(r)r J 

(vortex string) (13) 

and ks is the Boltzmann constant. We explain the notations in Eqs. (|1(])-(|13|) in terms 
of Coulomb-gas model. e(r) is the dielectric constant of the medium between unit charges 
separated by the distance r, and P(r) is associated with the density per area of thermally 
excited pairs of positive and negative charges (vortex and antivortex pairs). In Eq. fllPf), the 
increments of e(r) is given by the number of dipoles in the region surrounded by circles with 
radius r and r + dr. The term 2/(e(r)r) in Eqs. (|T2|) and ( |T3"D is the dipolar interaction term 
and the second term in Eq. (|l|) is the interlayer Josephson coupling. We have considered 
the screening only on dipolar interactions but not on the interlayer Josephson coupling. 
Dielectric constant e(r) is defined from renormalized coupling constant in the ab plane K{r) 
as 7re(r) = 1/K(r). The coupling constant K(r) is defined by 

K(r) = K\(r), (independent vortex) (14) 
K(r) = nKi(r). (vortex string) (15) 

The renormalized coupling constant for a single layer K\ (r, T) is defined as 

h 2 n 2D 

Ki(r,T)= , f (16) 

This equation is deduced from intervortex interaction energy as calculated by GL formalism 
in 2D. n 2 s D is the superconducting carrier density per area in a single 2D plane defined in 
terms of 3D carrier density rig as n 2D = n^d, m* is the effective mass of superconducting 
electron. The factor n in Eq. (p~5|) arises since the carrier density associated with vortex 
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string is n times that of independent vortex in a single layer. x(r, T) is the relative dielectric 
permittivity of medium in the 2D Coulomb-gas model in a single layer. For r = £ a &(T), 
we have x( r — £ab(T), T)=l and then x(r, T) represents the effect of screening by the 
renormalization described by Eqs. ( |T0| ) and (pTT|) . 

At r = r„, the differential forms of Eqs. (|It]) and ( |XT| ) are replaced by the continuity 
equations for e(r, T) and P(r, T). At first, we assume that x(r, T) is continuous at r = r n as 

X (r n + 0,T) = X (r n -0,T), (17) 

and also assume the conservation law for superconducting carrier density per area per layer 

n 2 s D {r n + 0, T) = nf (r n — 0, T) = const. (18) 

Here ±0 denotes the infinitesimal increment or decrement to r n . Combining Eq. (|18D with 
Eqs. (|1D-(|16|) yields the continuity of #i(r,T) as Ki(r„ + 0,T) = Ki(r n - 0,T) and the 
connectivity of K(r,T) as K(r n + 0,T) = nK(r n — 0,T). Thus, we have the connecting 
equation of e(r) as 

e(r n + 0,T) =e(r n -0,T)/n. (19) 

The connectivity of P(r) at r = r n is obtained as follows. A single independent vortex in 
a single layer is transmuted into a vortex string threading the film. Therefore, in n-layered 
film, the density of the vortex string per area in a film must be n times the density of 
independent vortex per area in a single layer. Thus, we have 

P(r n + 0,T) = nP(r n -0,T). (20) 

Putting the relation P(r n + 0,T) oc exp(— U st {r n + 0,T)/(fcsT)) and P{r n — 0,T) oc 
exp(— U vp {r n — 0,T)/(kBT)) into Eq. (PU[), we have 

U st (r n + 0,T) = U vp (r n - 0,T) - k B T\n{n). (21) 

The right-hand side of Eq. fl2~T| ) is F vp (r) of Eq. (fj). The Inn term of Eq. (^T|) arises from 
the entropy difference which in turn arises from the transition from independent vortex to 
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vortex string. These results suggest that we must compare the free energy of each excitation 
in selecting the favored excitations of n-layered system. 

Now, to solve the RG Eqs. ( |T0| ) and ([TT]) , we introduce the variables [^,[| 

' = ln fe) • (22) 

y%T)=r 1 P(r). (24) 

K(l,T) is used in place of K(r,T). By introducing the variables I, x(l,T), and y(l,T), we 
can rewrite Eqs. (4) and (5) in terms of Coulomb-gas model in a picture of vortex. Thus, 
the RG equations (4) and (5) are transformed into the form, 

dx{l > T) =8**y%T), (25) 



dl 

dy(l,T) 2y(l,T) 



dl x(l,T) + l 



x 



H- 2 ^^ (2i) ). 

2D KT RG equations of Eqs. fl2"5| ) and are written with respect to independent vortex 
and antivortex pairs in a single layer for I < l n and with respect to vortex-string and 
antivortex-string pairs threading the film for I > l n . These equations are written in the 
same form except the coupling constant K(l,T) and K±(l,T) defined as follows, 

, KS,T), (I < l n ) 
K(l,T) = { U ' h 1 n) (27) 

nK x {l,T), (l>l n ) 
, + 0, (1 < Q 

KdhT) r ( 28 ) 
= 0. (l>l n ) 



The second term in the right-hand side of Eq. (|26| ) for I < l n is an additional term to the 
usual 2D KT RG equations and represents the interlayer coupling. Introducing the variable 
for a single layer Xi(l,T) = 7r ^-^ 7 r - ) — 1, we have the relation 
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x{l,T)= Xl {l,T), (/</„) (29) 
x(l, T) + 1 = (xi(Z, T) + l)/n. (/ > Z„) (30) 

The initial conditions for x(0,T) and y(0,T) in the RG differential equations (p5|) and 
(p6|) are derived from the fact that initially at Z=0 only independent vortex and antivortex 
pairs in each layer exist for the case K±/K\ <Cl. The initial condition will be reconsidered 
in detail in Sec. 5 taking account of the entropy corrections in the system. Then, x(0,T) is 
given by |,| 

«W=M = (3D 

Here, Ki(0,T) is written in terms of magnetic penetration length X a b{T) as K\(0, T) = 
(h 2 dc 2 eo)/(Ae 2 Xl b (T)). |7||| eo is the dielectric constant of vacuum, e is the unit charge, and 
c is the velocity of light. We assume X a b(T) = A o &(0)(l — T/T c0 (n)) _ 5 ? and A a b(0)=260 nm. 
|T8| Now, we rewrite Eq. fl3T| ) as 

The initial value for y(0,T) is given in terms of the initial value x(0,T) and the parameter 
B by 

2x(0, T) - 2 In (x(0, T) + 1) — 4ttV(0, T) 

= -fi(T/T^ T -l). (33) 

Following 2D KT theory, we assume the linear T dependence in the right-hand side of 
Eq. fl33|). Here, the parameter S is defined as [| 



1 T* 

5 ~ - (34) 

6T c0 (n)-T£ T ' 1 ; 

In our model, the parameters b and T c0 (n) are determined from the experimental data as 
listed in Table |. We have used single layer T^ T for KT-transition temperature in Eq. ([3T]) . 
Tk T is taken as a scaling parameter for T to determine the initial values x(0, T) and y(0, T) 
for a given temperature T. 
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The connection between the differential equation of the two regions, i.e., the region I < l n 
and the region I > l n , is expressed as continuities of Kx(l,T) and Xi(l,T) at I = l n . Then, 
the connectivity condition in terms of x(l,T) is written as 

x(l n + 0, T) + 1 = (x(l n — 0, T) + l)/n. (35) 

The connectivity condition for y(l,T) is written as 

y(l n + 0,T) = ,/Ey(l n -0,T). (36) 



Equation fl36|) is obtained by inserting the connectivity of P(r, T) of Eq. fl20|) into Eq. 

At the KT-transition temperature, and for the limit I — > oo, we have x(l,T' c (n)) — > and 

y{l,T' c {n)) — > 0. Then, from Eq. (|30D, we have xx(l,T' c (n)) — > n — 1. 

Using the parameters of BSCCO films listed in Table p], the RG differential equations of 
(p5|) and (^) for x(l,T) and y(l,T) were solved with respect to the variable Z for the case 
n=l, 3, 6, and 20. The parameters l n {T' c {n)) are derived, using the parameter b, T c0 (n), 
and renormalized K\/K± at T' c (n). The parameters T c o(n), and 6 are determined from 
the experiment as will be discussed in Sec. 4 (Fig. 0). The initial value of the anisotropy 
parameter Kx/K± is obtained as K\j K± = 10 3 using the parameters used in Table |. 
However, it is reduced from the initial value by renormalization as I increases. Renormalized 
Kx(l, T) is used in each step of discretized I during the numerical calculation. We will confirm 
in Sec. 5 that the lowest free-energy excitations are such as we have assumed initially even 
after the renormalization of Ki(l,T). 

The solutions of the RG equations are shown in Fig. |l| The RG flow for the case n=l is 
shown in Fig. |l|(a). The parameters b and T c o(l) are extrapolated to the case n—1 from those 
for n=3, 6, and 20 shown in Table |. The points (x(0, T), 2-7n/(0, T)) are shown as a dashed 
curve. The points (x(l,T),2iry(l,T)) follow the flow of 2D XY model as I increases. In our 
numerical solution the RG equation of Eq. ( |26| ) includes the factor l/(x(l,T) + 1) on the 
right-hand side of Eq. (|2"6]). In the limit | x(l, T) |— ► 0, the solution follows a usual parabolic 
one. The RG flow starts from a;(0,T^(l)) for T c '(l)/T£ r =1.0, and the point converges 
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to x(l,T'J\)) -> 0, and 2<Ky{l,T' c {l)) -> 0. For T > T c '(l), 2iry(l,T) diverges hyperbolically 
and then asymptotically linearly to infinity in the limit x(l,T) — > oo. For T < T c '(l), x(l,T) 
is confined within limited range for 2ny(l,T) > 0. 

RG flows for multilayered thin films are shown in Figs. 1(b)- (d). We select the sin- 
gle layer KT-transition temperature as T^ T =30 K. In these figures, we show the contour 
of points (xi(l, T), 2Try(l, T)) for I < l n and (xi(l,T), 27n/(Z,T)/y / n) for I > l n instead of 
(x(l,T),2ny(l,T)) in the case of n—1. Initially for I <C l n the points (x\(l, T), 2iry(l, T)) 
are associated with the pairs of independent vortex, and follow the flow of pure 2D XY 
model. As I increases further, the RG flow is deflected downward, i.e., density of indepen- 
dent vortex is largely reduced, owing to the influence of the r 2 interlayer coupling term. 
Then, at I = l n the flow shows a transition to the region of vortex string. The transi- 
tion occurs continuously or shows a slight cusp depending upon temperature. The points 
— 0,T), 2vn/(Z n (T) - 0, T)) or equivalent^ {x l {l n (T) + 0, T), 2vn/(/ n (T) + 0, T)/y/n) 
for various temperatures are connected as dashed curves in Figs. l(b)-(d). When I > l n , 
the points (xi(l, T), 2Try(l, T)/y/n) follow the 2D XY RG flow with respect to vortex string. 
Below and above the temperature T' c (n), the RG flow behaves differently. In the region 
T < T' c {n), the flow has limited range, while for T > T'Jji) the flow diverges parabolically 
and then asymptotically linearly like 2iry/y/n = x/y/n to infinity following the 2D KT RG 
flow with respect to vortex string and the density of dissociated vortex string increases as 
a function of I. This is the KT transition with respect to vortex string under the influence 
of interlayer coupling. From this, we have the transition temperature T^(n)/T£ T =1.47±0.1, 
1.58±0.1, and 2.16±0.1 for the cases n=3, 6, and 20, respectively as listed in Table |. Note 
that xi(l,T) converges to n — 1 as y(l,T) — > at T=T' c (n). The renormalization flow 
indicates that the vortex string persists below and above T' c (n). 
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III. SAMPLE PREPARATION AND MEASUREMENT 



Samples of BSCCO thin films with various thicknesses were prepared by MBE on 



Nd:YA103 (001) substrates, as described elsewhere. |19-21| All of the films are c-axis ori- 
ented. Current flow within the ab plane was measured so that the circular current, which 
forms the vortex aligned along the c axis, are the cause of the resistive transition. Usually, 
BSCCO thin films are twinned with respect to the a and b axes. On the contrary, our films 
grown on Nd:YA103 (001) are untwinned. |T9"^PT| We used the untwinned films to reduce 
the pinning effect of the vortices by the twin boundaries. We have obtained the exact film 
thickness by the following relation; 1.54 nm x number of stacked layers. The stacked number 
is estimated from the oscillations due to the Laue function in the x-ray diffraction pattern. 
These samples showed good superconducting transitions. 

Electrical measurements were done by D.C. current method with a usual four-terminal 
method. The four contacts ( Au wire with diameter of 25 /im, Ag paste) were made on the 
ab plane of the sample. During the measurement, care was taken to keep the temperature 
constant within the error of 0.05 K. Errors in the indicated temperatures were kept within 
0.1 K by minimizing the difference in measured electrical characteristics obtained during 
cooling and heating cycles. 



IV. EXPERIMENTAL RESULTS AND KT TRANSITION BY VORTEX STRING 

In this section, temperature dependence of resistance R and nonlinear I — V characteris- 
tics on thin films at various temperatures near T c are described. These results are interpreted 
in terms of KT scaling by vortex strings. Films with their thicknesses less than ro should be 
used in order to eliminate vortex rings. To obtain the experimental parameters to resolve 
the RG equations, we used the data of the film with n=20, which is thin enough to exclude 
the effect of the vortex ring. The typical resistance behavior is shown in Fig. 2. The sam- 
ple size along the ab plane is 2.0 mm x 6.0 mm and the thickness is 30.9 nm. The inset 
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of Fig. 2 shows a representative result of the electrical resistance of thin films with n=20 
as a function of temperature. This figure shows the superconducting transition tempera- 
ture (zero-resistance temperature) T c (20)=64.5 K, the onset temperature T onset ~79 K, and 
the midpoint temperature (20) =71. 5 K, which is used as the GL mean-field transition 
temperature T c0 (20). 

Applying the pure 2D KT theory to vortex string, we have the resistance R that 

arises from the flow of dissociated free vortex strings at T > T' c {n) as 



T cQ (n) - T 



R/Rn oc exp ^— 2 A 

Here, Rn is defined by the normal-state resistance at the onset temperature. Note that T' c (n) 



is used in place of KT-transition temperature for pure 2D KT theory. Equation (|37| ) is based 
on the assumption that KT transition occurs owing to the dissociation of vortex-string and 
antivortex-string pairs as we have shown in Sec. 2. 

Normalized resistance R/Rn is shown on a logarithmic scale as a function of the nor- 



malized temperature y ^Zj^jy i n Fig- 2. Data points lie on a straight line around T c (20). 
The coefficient b = 1.7 is derived from the gradient of the R vs T curve, using T c (20)=64.5 
K, and T c0 (20)=71.5 K in Eq. (j37|). The n dependence of b is rather small as is shown in 
Table |. We extrapolated these values of b for n > 1 to the value 6=1.7 for n=l in olving 
RG equation as was described in Sec. 2. 

According to KT theory, I — V characteristics should exhibit nonlinear power-law be- 
havior V oc J a ( T ) owing to the current-induced dissociation of the pairs of vortex string and 

antivortex string. The exponent a(T) of vortex string is given by [0,^,0 

ttK(LT) ixtfn 2 ? , , 

a T = 1 + l ' J = 1 + ' 38 

KbI 4m*x{r, 1 )Kb-l 

This is based upon our RG analysis where the vortex string persists below and above T' c {n). 
According to KT theory, the exponent oi(T) — 1 jumps from 2 to at T' c (n). This point 
determines T' c {n) as we shall discuss in detail in Sec. 5 (See Eq. (ffSD). 

To see the variation of oc(T), the I — V characteristic was measured at various temper- 
atures around T c (20), using the same sample, as is shown in Fig. 3. The I — V curves on 
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a log- log plot yield straight lines over 2 orders of magnitudes in current, indicating power 
law V oc I a<yT \ The gradients of the lines get steeper as the temperature gets lower around 
T c (20). Figure 4 shows the exponent oi(T) — 1 as a function of temperature obtained from 
the I — V characteristics in Fig. 3 and the renormalized coupling constant (tvK(1,T)) /(ksT) 
(namely l/(x(l,T) + 1)) calculated in Sec. 2. The filled circles show the measured points. 
These circles are compared with the solid curves of the calculated coupling constant. Be- 
low 66 K, the exponent a(T) — 1 jumps from to 2 and increases more steeply at lower 
temperatures. At T ~ T c (20)=64.5 K, the exponent a(T) — 1 ~ 2 holds. 

The coupling constant (ttK (l,T)) / (ksT) of Eq. (|38|) is drawn as a function of T around 
T c '(20) for 1=3, 4, 5, 10, and 20 in Fig. 4. These / regions are, as we have mentioned in 
Sec. 2 (Table [I]), in the vortex-string region. For a given I and the temperatures higher 
than T c '(20), (ttK(1, T))/{k B T) drops abruptly. For the larger values of I, the gradient gets 
steeper. In the T range below and near T c '(20), the theoretical curves can be fitted to the 
experimental results within an experimental error. Above T c '(20), the calculated coupling 
constant for the vortex strings of 1=4 fits the experimental curve best. Then, we can estimate 
the size of interacting vortex-string pairs as r = £ a f,(T) exp (3) ~590 nm, near the transition 
temperature T c '(20). 

The various transition temperatures T c (n), T' c {ri), and T c o(n) are shown as functions of 
n in Fig. 5. To calculate T c '(ra), we have used the experimental values for T c0 (n), b, X a b(T), 
as are listed in Table |, and K 1 /K± = 10 3 (initial). Here T mid (n) is regarded as T c0 (n). The 
agreement between the experimental T c (n) and the calculated T c '(n) are satisfactory. 

V. DISCUSSION AND CONCLUSION 

Using the solution of the RG equations, the energies of U vp (r), U r i ng (r), U st {r) of Eqs. ([[])- 
(0) and the free energies of F vp (r) of Eq. (^) are shown as a function of r for the case n=20, at 
temperatures around T c '(20) in Fig. 6. The dashed curve representing F vp (r) intersects the 
curve of U st {r) at r = r n . We used the renormalized Ki/K± obtained in Sec. 2 in Eqs. ([!])- 
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(0). This figure clearly visualizes that the lowest free-energy excitation changes from that 
of independent vortex pair to vortex-string pair at r = r n at all temperatures. At large r, 
the free energy of a vortex-string pair is the lowest because of its weakest r dependence. 
At T = T c '(20), the free energy of a vortex-string pair increases slowly like lnr because the 
coupling constant becomes a constant value K\ ~ 2k B T' c {2Q) /-n. At T > T c '(20), it declines 
because K\ decreases rapidly. At T < T c '(20), although not shown, it behaves similarly as 
it does at T = T c '(20). For the characteristic length r n the renormalized Ki/K± is used. 
Notice that the relation nd < r n still holds after renormalization. A cusp of the energy 
of vortex-string pair at r = r n appearing in Fig. 6(b) for T > T^(20) is due to the large 
change of Ki(l,T) at the transition from independent vortex pairs to vortex-string pairs 
that appeared in RG flow in Fig. 1(d). The position of r=590 nm, at which the a{T) — 1 vs 
T relation fits the experimental result best, is located outside the range of Fig. 6. 

Now let us reconsider the initial condition for RG equations of Eqs. (^) and fl2"6|) in 
some detail by taking account of entropy corrections. We need some justification for this 
condition assumed in Sec 2. We are discussing a narrow region near r « £, a b{T) in Fig. 6 
where energies of the three excitations are close to one another. We have a narrow region near 
r= £,ab{T), although not recognizable in Fig. 6, where vortex rings are the lowest free-energy 
excitations but independent vortex pairs are not for the case £, a b{T) < nd and K\/K±_ > K r . 
For the case £, a b(T) > nd and Ki/Kj_ > K p , we have a narrow region near r=£ a (,(T) where 
vortex-string pairs are the lowest free-energy excitations but independent vortex pairs are 
not. However, for large anisotropic ratio, the widths of these regions are of the order of 
U{T){KjK x ){U{T)/df or U{T)^{KJK{) or £ a b(T) (Kx/ Ki) (£ab(T)/ d) . In addition, 
the energy difference among excitations is also small in these regions. So if we regard that 
the initial lowest free-energy states are independent vortex for solving RG equations, errors 
should be small. 

Now let us discuss T' c {n) vs n relation obtained in Sec. 2. T' c (n) is defined by the 
temperature where x(l,T{,(n)) — > and 7/(/,T c '(n)) — > are satisfied as / — ► oo. Accordingly, 
from Eq. (p3|) and using the relation Ki(0, T' c (n)) oc (1 — T^.(n)/T c0 (n)), we have 
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where the relative dielectric permittivity x ( oo, T c '(n)) is defined as T' c {n)) = 

Ki(Q, T c / (n))/iT 1 (oo, T' c {n)). T\ is the scaling temperature for T c0 (n) defined as T\ = 
(nh 2 dc 2 e )/(8k B e 2 \ 2 ab (0)). For the case of n=l, Eq. © yields T' c (l)=T* KT . Using the 



parameters for BSCCO ||], we have T A =223.7 K. From Eq. ©, T c0 (n) and T£ T set the 
upper and lower bounds for the T' c {n) vs n relation, respectively. The experimental result of 
T c (n) vs n relation shown in Fig. 5 is similar to that of YBCO in a sense that T c0 (n) — T^ T 
is of the same order of T^ T itself. This means that the interplane interaction is equally 
important as the intraplane KT interaction. The difference lies in the shape of T' c (n) vs n 
relation. This difference arises from the difference in n dependence of x(°°> T' c {nj), which in 
turn originates from the difference in the anisotropy ratio K\ / K± , as is indicated in Eq. ( ft^D . 
For smaller anisotropy parameters, i.e., Ki/Kj_ «1, l n gets smaller, and intervortex inter- 
actions are screened largely by vortex string pairs. In this case, T' c {n) approaches that of 
2D KT theory with respect to vortex string. Matsuda et a/.'s model corresponds to the case 
of the limit l n =0 in our model. In contrast, for larger anisotropy ratio K\jK\_ 3>1, l n gets 
larger and intervortex interactions are screened by independent vortex pairs with r < r n 
and by vortex-string pairs with r > r n . In this case, T' c {n) approaches less rapidly that of 
2D KT model with respect to vortex string. This is the case for BSCCO and to solve the 
RG equation is inevitable to obtain T' c {n) vs n relation. 

Now, within the framework of our model, we discuss the absence of dissociation for 
independent vortex and antivortex pairs. In our model, KT transition occurrs in vortex- 
string region. This is due to the absence of r 2 term as is indicated in Eq. (£J) in this region. 
This can be understood by the following reasoning. If we have the interlayer coupling term 
in the form of K±(R C / 'dy '(7 > 0), the free energy of the system with an isolated independent 
vortex is given by F SV (R C ) = nK\ ln(R c /£) + K±(R c /d)' y — 2kBT\ia(R c /^). The last term is 
the entropy term and R c is the sample size along the ab plane. In the limit R c — > 00, we 
have F SV (R C ) — > 00, irrespective of T. This means the absence of vortex dissociation. This 
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is also suggested by the RG flow in Fig. 1. The RG curves seem to continue to fall off at 
large I region if transitions to vortex string were absent. Thus, in general, we have bound 
state for infinite sample as far as interlayer coupling is a more steeply rising function of r 
than lnr. Therefore, the low-lying excitation should necessarily be other than independent 
vortex pairs in infinite area thin films, i.e., vortex-string and antivortex-string pairs in thin 
films. 

In the case of finite and very anisotropic samples with limited area along the ab plane, 
however, KT transition by independent vortex and antivortex pairs may occur. Equating 



the first and the second term of Eq. (|l|), we have r& = dJ(Kx/K±)]xi(rb/^ a b(T)). [0 For 

r < r&, the first lnr term is larger than the second r 2 term. The length becomes large 
as the anisotropy ratio Kx/ K± becomes large. If the sample size along the ab plane is less 
than rfe, r 2 -dependent coupling term is negligible as compared with lnr term within the r 
range of the sample size and 2D KT transition can occur. In this case, however, T c (n) will 
be independent of n but depend on sample size along the ab plane when this size is near 
rfe. We believe that this case is irrelevant for our very thin films of BSCCO, since T c (n) 
increases as a function of n. 

Finally, we will make some comments on the significance of our consequence described so 
far in bulk high-T c superconductors. Following Matsuda et al, we assign the superconducting 
transition temperature of bulk BSCCO crystal to 2D/3D crossover temperature T*. This 
is defined by Hikami and Tsuneto as the temperature, at which x(r ~ ?"f>,T*) diverges and 
given by f| 



In RG flows such as in Fig. 1, the position of x(r ~ r^T) where the RG flow starts to fall 
off goes to infinity at T = T*. In our RG flows of Fig. 1 with temperatures around T' c (n) 
these deflecting points are located within finite regions. This means that the calculated 
superconducting transition temperatures T' c (n) are less than T*. These results are consistent 
with our view that the T c of our thin film is due to vortex string and that of bulk is 2D/3D 





(40) 
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crossover. If we insert T* ~70 K and T£ T =30 K, we have (T* - T^ T )/T^ T ~1.3. As far as 
we understand, we should use the renormalized K\j K± in Eq. ( f40|) as we have obtained l n 
by solving RG equations. For example, at T=70 K, we have 7"b=28 nm, Ki(rb) / K±=160, 
and the right-hand side of Eq. ( f40|) reduces to 1.5. A small discrepancy may be due to the 
renormalization of K± in the bulk. 

On experimental side, we have grown very thin films by MBE and verified vortex-string 
excitations therein. We have used untwinned films exhibiting sharp superconducting transi- 
tions as described in Sec. 3. However, more systematic experiment determining the physical 
parameters of thin films are needed to quantify the model. In particular, superlattices with 
superconducting layers sandwiched by nonsuper conducting layers are very interesting to 
study the role of interlayer coupling. The measurement of T c (n) vs n relation on other high- 
T c thin films such as L^-zSr^CuO^ and TBCCO is important for the firm establishment of 
the concept of vortex string. 

So far, we have interpreted our experimental results on superconducting properties of 
BSCCO thin film in terms of KT transition by vortex string. In this view, the influence of 
pinning centers such as defects can be observed as a perturbation to the KT transition. The 
main reasons that we interpret our experiment by KT transition are following. (1) There is 
a broad transition region in superconducting transition. The transition temperature T c is 
far lower than T c0 . (2) Universal jump from 3 to 1 is observed in a-T relation. (3) The n 
dependence of T c is observed in films of various numbers of layers. The T c -n relation is the 
experimental evidence of vortex string. In addition, it is believed that the pinning effects 
are weak in BSCCO superconductors ]22| and that our BSCCO film grown on Nd:YA103 



(001) substrates are untwinned as was described in Sec. 3. If the transport properties of 
vortices are not dominated by vortex-antivortex interaction but by vortex-defect interaction, 
these behaviors characteristic of KT transition will not be observed. For example, in the 



presence of high-density defects, the universal jump will be obscured. |23[ We showed by our 
RG scheme that to include r 2 dependent interaction necessarily leads to the KT transition 
by vortex string. Comparing the solution of RG equation with the experimental results 
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of a-T relation, we obtained I— 4, the logarithm of the size of interacting pairs in the unit 
of £ab(T) as was described in Sec. 4. In our view, the finite size of independent vortex 
and antivortex pair or vortex-string and antivortex-string pair represents the distance in 
which independent vortex and antivortex or vortex-string and antivortex-string can migrate 
without being influenced by the defects. From the a-T relation for samples with n=3, 
we have / ~2.5, shorter than 1=4 for n=20. This means that the mean free pass for vortex 
transport becomes shorter as the film becomes thinner and T c (n), lower. In this case, pinning 
effects are stronger and vortices are more firmly pinned. This tendency seems physically 
reasonable, since activation energies for vortices to be freed from pinning are lowered at 
higher temperatures. Note that the value / ~2.5 is rather close to Z n =2.34, which suggests 
that pinning effect near T c (n) is fairly strong in the films n=3. In this way within the 
framework of our view, the effects of the pinning can be considered as a perturbation to the 
KT transition by vortex string and assessed consistently in terms of finite length of I. Thus, 
we believe that our view for our experimental results is basically correct. 

In conclusion, we have formulated the KT theory with interlayer coupling to quantify 
the concept of vortex string. We write down the KT RG equations on independent vortex 
pairs for r < r n and on vortex-string pairs for r > r n and then set up the scheme connecting 
the two regions. KT transition occurred in the vortex-string region with lnr interstring 
interaction. We believe that the fundamental assumptions in our model are basically correct 
since I — V characteristics and the T' c (n) vs n relation in our BSCCO films were interpreted 
by vortex string. Significance of these results was discussed in brief in relation to bulk 
high-T c superconductors. 
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FIGURES 

FIG. 1. Renormalization contours of vortex and antivortex pairs. For the case (a) the points 
(x(l,T),2iry(l,T)) are shown. For the cases (b)-(d), the points (xi(l,T),2iry(l,T)) in the range 
/ < l n are shown and (xi(l,T),27ry(l,T)/yfn) are in the range I > l n . The points correspond to 
intervals of Z=0.1. The parameter 6=1.7 is used. For (b)-(d) Ki/K±=W 3 is used, (a) n=l, 
and K±=0. The curves are with T/T£ T =0.5, 1.0, and 1.3 (from left to right). The initial 
points (x(0, T), 27ry(0, T)) are shown as a dashed curve, (b) n=3. T/T^ T =1.40, 1.47, and 1.49. 
(c) n=6. T/T£ r =1.55, 1.58, and 1.62. (d) n=20. T/T£ T =2.13, 2.16, and 2.20. The transi- 
tion points from independent vortex pairs {x\{l n — 0,T),2iry(l n — 0,T)) to vortex-string pairs 
(xi(l n + 0,T),27ry(l n + Q,T)/y/n) for n=3, 6, and 20 (Fig. |l|(b)-(d)) are connected by dashed 
curves. 

FIG. 2. Normalized resistance R/Rn as a function of normalized temperature. T c (20)=64.5 
K and T c o(20)=71.5 K are obtained from the inset. Rn is defined by the normal-state resistance at 
the onset temperature. Inset: Temperature dependence of the resistance for BSCCO film (n=20) 
at / = 3.5 x 10" 6 A. 

FIG. 3. Current-voltage characteristics at various temperatures for BSCCO film (n=20). 
Numbers in the figure denote temperatures. 

FIG. 4. The a(T) — 1 obtained from the I — V characteristics as a function of temperatures in 
Fig. 3 and the calculated curves of the renormalized coupling constant {nK{l,T)) / (ksT) (namely 
1 /(x(l, T) + 1)) derived from Fig. |l|(d). The filled circles show the T dependence of measured points 
of the exponent a(T) — 1. The solid curves are theoretically calculated renormalized coupling 
constant. 
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FIG. 5. Experimental data of T c (n) and T c q(tl), and the calculated values T' c (n) as a function 
of the number of layers n. Filled circles are the experimental data of superconducting transition 
temperature T c (n), open circles are the experimental data of GL mean-field transition temperature 
T c o(ra), when we assume T c o(n) = T m id{n). The open squares are the superconducting transition 



temperature T' c (n), which are obtained from the RG equations (25) and (^3). Data points of T c (n) 
and T c g(n) for n=3 and 6 half unit-cell thin films are taken from Refs. 20 and 21. Solid, dot-dashed, 
and dashed curves are guides to the eyes for T' c (n), T c (n), and T c o(n), respectively. 

FIG. 6. The energies of three types of excitations, U vp (r), U r i ng (r), and U s t(r) and the free 
energy F vp {r) as a function of r in the unit of k B T for n=20. (a) T = T c '(20). (b) T > 1^(20) 
(T/T^ T =2.20). U(r) are shown as solid curves and F vp (r), as dashed curves. The renormalized 
K\(r,T) is used. Their intersections at nd and r n are shown. 
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TABLES 

TABLE I. List of parameters for BSCCO thin films. £a&(T c (n)), l n (T' c {n)), b, T c o(n), T m id(n), 
T c (n), T' c (ri), and T' c {n) /T^ T for n=l, 3, 6, and 20. £ a £>(T c (n)) are calculated using £ a £>(0)=3 nm and 
experimental values of T c (n) and T c o(n). l n {T' c (n)) are calculated using Eq. (|2l| ) and renormalized 
Ki(Z n (T),T), and d=1.54 nm. The initial value of K x (l = 0,T)/K±(l = 0,T)=10 3 is calculated 



using £ c (0)=0.1 nm. [24] b for n=20 is taken from R vs T relations in Fig. 2. T c o(n) or T m id(n) 
for ra=20 is taken from the inset of Fig. ||. T c (n) for n=20 is determined from R vs T relations in 
the inset of Fig. T c '(n) is calculated from RG equations (^) and ( |2"6| ) using experimental values 
listed in this Table. T£ T is taken as T^- T =30 K. The parameters T c (n), T c o(n), and b for n=3 and 
6 are taken from Refs. 20 and 21. T c o(n) and 6 for n=l are extrapolated from those for n=3, 6, 
and 20. For the definition of parameters, see the text. 

n 1 3 6 20 

U(T c (n)) (nm) 4.2 4.8 5.2 9.6 

l n {T' c {n)) (0) 2.34 2.54 2.70 

b (1.7) 1.4 1.7 1.7 

T c0 (n),T mid (n) (K) (62) 62 65 71.5 

T c {n) (K) 38 43 64.5 

T' c (n) (K) (30) 44.1±3 47.4±3 64.8±3 

T' c {n)lT* KT (1) 1.47±0.1 1.58±0.1 2.16±0.1 
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